ON COMPUTABILITY OF JULIA SETS: ANSWERS TO QUESTIONS OF 

MILNOR AND SHUB 

M. BRAVERMAN AND M. YAMPOLSKY 

\o ■ 

^— ^ , Abstract. In this note we give answers to questions posed to us by J. Milnor and M. Shub, 

f^ ' which shed further light on the structure of non-computable Julia sets. 
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1. Introduction 



Computability of real sets. The reader is directed to |BYj for a more detailed discussion of 
the notion of computability of subsets of M"" as applied, in particular, to Julia sets. We recall 
the principal definitions here. The exposition below uses the concept of a Turing Machine. 

^ \ This is a standard model for a computer program employed by computer scientists. Readers 

^ • unfamiliar with this concept should think instead of an algorithm written in their favorite 

programming language. These concepts are known to be equivalent. 

Denote by D the set of the dyadic rationals, that is, rationals of the form ^. We say 
that : N ^ D is an oracle for a real number x, if \x — (f){n)\ < 2~"' for all n & N. In 
other words, provides a good dyadic approximation for x. We say that a Turing Machine 

C^- ■ (further abbreviated as TM) M'^ is an oracle machine, if at every step of the computation 

M is allowed to query the value (f){n) for any n. This definition allows us to define the 

(3 ' computability of real functions on compact sets. 

^' 

O I Definition 1.1. We say that a function / : [a,b] -^ [c,d] is computable, if there exists an 

^ ■ oracle TM M'^{m) such that if is an oracle for x G [a, 6], then on input m, M'^ outputs a 

t^ : 1/ e D such that \y - f{x)\ < 2""^. 

It is worthwhile to note why the oracle mechanism is introduced. There are only countably 
many possible algorithms, and consequently only countably many computable real numbers 
/\ ' which such algorithms can encode. Therefore, one wants to separate the hardness of encoding 

j^ ■ the real number x from the hardness of computing the value of the function /(x), having 

the access to the value of x. 

Let iiT C M*^ be a compact set. We say that a TM M computes the set K if it approximates 
K in the Hausdorff metric. Recall that the Hausdorff metric is a metric on compact subsets 
of M" defined by 

(1) dffiX, Y) = inf{e > \ X C U,{Y) and Y C U,{X)}, 

where Ue{S) is defined as the union of the set of e-balls with centers in S. 

We introduce a class C of sets which is dense in metric dn among the compact sets and 

which has a natural correspondence to binary strings. Namely C is the set of finite unions 
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of dyadic balls: 

r " 

C = < [jB{di,ri) I where di,ri G 1 

Members of C can be encoded as binary strings in a natural way. 

We now define the notion of computability of subsets of M" (see |Weij . and also |RWj ) . 

Definition 1.2. We say that a compact set i^ C M'^ is computable, if there exists a TM 
M{d,n), where rf G ©, n G N which outputs a value 1 if dist{d,K) < 2"", the value if 
dist((i, K) > 2 ■ 2"", and in the "in-between" case it halts and outputs either or 1. 

In other words, it computes, in the classical sense, a function from the family J^k of 
functions of the form 

if dist{d,K) >2-2-" 
(2) f{d,n) = { I, if dist{d,K) < 2-" 

otherwise 

Theorem 1.1. For a compact K cM.'^ the following are equivalent: 

(1) K is computable as per definition M/A 

(2) there exists a TM M{m), such that on input m, M{m) outputs an encoding of Cm G C 
such that dniK, Cm) < 2^*" (global computability), 

(3) the distance function dxix) = mi{\x — y\ \ y G K} is computable as per definition \l.l\ 

Note that in the case k = 2 computability means that K can be drawn on a computer 
screen with arbitrarily good precision (if we imagine the screen as a lattice of pixels). 

In the present paper we are interested in questions concerning the computability of the 
Julia set Jc = J{fc) = J^z^ + c). Since there are uncountably many possible parameter 
values for c, we cannot expect for each c to have a machine M such that M computes Jc 
(recall that there are countably many TMs). On the other hand, it is reasonable to want M 
to compute Jc with an oracle access to c. Define the function J : C — »• K* [K* is the set of 
all compact subsets of C) by J{c) = J{fc)- In a complete analogy to Definition II. II we can 
define 

Definition 1.3. We say that a function n : S ^ K* for some bounded set S is computable, 
if there exits an oracle TM M'^{d,n), where (j) is an oracle for x E S, which computes a 
function (0) of the family J-'k(x)- 

Equivalently, there exists an oracle TM M'^im) with again representing x E S such that 
on input m, M''^ outputs a C G C such that dfjiC, k{x)) < 2~™. 

In the case of Julia sets: 

Definition 1.4. We say that Jc is computable if the function J : d \-^ J^ is computable on 
the set {a}. 



We have the following (see |BBY2j ): 



Theorem 1.2. Suppose that a TM M"^ computes the function J on a set S G C Then J is 
continuous on S in Hausdorff sense. 
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Previous results. We have demonstrated in IBYI: 



Theorem 1.3. There exists a parameter value c G C such that the Julia set of the quadratic 
polynomial fc{z) = z^ + c is not computable. 

The quadratic polynomials in Theorem 11.31 possess Siegel disks. It was further shown by 
I. Binder and the authors of the present paper in |BBYlj that the absence of rotation 
domains, that is either Siegel disks or Herman rings, guarantees computability of the rational 
Julia set. This implies, in particular, that all Cremer quadratic Julia sets are computable - 
this despite the fact that no informative high resolution images of such sets have ever been 
produced. One expects, however, that such "bad" but still computable examples have high 
algorithmic complexity, which makes the computational cost of producing such a picture 
prohibitively high. 

Two questions on computability of Julia sets. J. Milnor has asked us the following 
natural question: 

Is the filled Julia set of a quadratic polynomial always computable? 

In this paper we answer in the affirmative: 

Theorem 1.4. For any polynomial p{z) there is an oracle Turing Machine M^{n) that given 
an oracle access to the coefficients of p{z) and n, outputs a 2""^ -approximation of the filled 
Julia set Kp(^z)- 

Moreover, in the case whenp{z) = z'^+c is quadratic, only two machines suffice to compute 
all non-parabolic Julia sets: one for c G M, and one for c ^ M. 

This may come as a surprise, given the negative result for Julia sets. To gain some insight 
into how non-computability can be destroyed by filling in, consider the following toy example. 
Let y4 : N ^ {0, 1} be any uncomputable predicate. Consider the set 

^1 J {^e2-/fc|^e[i-i,i]} ioTt = {0.d,d2ds...)2e[0,i) 

Si \j {^e2^i/fc|^ e [l-i,l]} fort = l 

fcGN, A(fc)=l 

To avoid ambiguity, we always take the finite expansion for dyadic t's. An example of a set 
Qt is depicted on Figure H Firstly, note that if t G (0, 1) is not a computable real, then the 
set Qf is non-computable by a TM without an oracle for t. Moreover, even for a TM M'^' 
equipped with an oracle input for t, the set fli is clearly non-computable. However, when 
filled, every fit becomes a computable set - the unit disk. 

The question of M. Shub again has to do with fragility of non-computability. This time, 
instead of filling in a non-computable Julia set, we will make a "fuzzy" picture of it by 
letting the parameter c vary in a neighborhood. 

To formalize this, consider the following definition. Let J be the subset of C x C given by 



J= {{z,c) : z G Jc}. 
Shub has asked us: 
Is the set J computable? 
The answer again is "yes" : 



2m/2__i 




Figure 1. Part of the picture of Qt for t = (0.10101 . . .)2 
Theorem 1.5. Let d > be any computable real. Then the set 



jnCx B{0,d) 
is a computable subset ofCxC. 

Informally, we may think of projection of JI fl C x (c — e, c + e) to the first coordinate as the 
picture that a computer would produce when J^ itself is uncomputable. 

To understand how the mechanism of non-computability is destroyed in this case, consider 
again the set fit for t G (0, 1] as the toy model. The set W = {{z,t) : z e ^t, t ^ (0, 1]} C 
C X M is computable even though Qt itself is non-computable for t = 1. This happens because 
in the closure of W the "slice" corresponding to t = 1 is 

S^[j{re^^'/''\re [1 - ^, 1]} D f^i. 

rC 

fcGN 

This set "masks" the computational hardness of Qi, and makes W computable. 

2. COMPUTABILITY OF FILLED JULIA SETS 

Our goal in this section is to prove Theorem ll.4l For a given polynomial p{z) we construct 
a machine computing the corresponding filled Julia set Kp. We will use the following combi- 
natorial information about p in the construction. Note the this information can be encoded 
using a finite number of bits. 

• Information that would allow us to compute the non-repelling orbits of the polynomial 
with an arbitrary precision. Note that there are at most degp — 1 of them. Such 
information could, for example, consist of the list of periods ki of such orbits; and 
for each i a finite collection of dyadic balls {D^}j4i separating the points of the 
corresponding orbit from the other solutions of the equation p''^{z) = z. This allows 
for an arbitrarily precise approximation of the orbits by using an iterative root-finding 
algorithm for p^*{z) = z in Dj. 

• In the case of a hyperbolic or a parabolic orbit (, a domain of attraction D^ such 
that every orbit converging to ( eventually reaches and stays in D^. In the hyperbolic 
case D(^ is just a collection of discs. In the parabolic case, it is a collection of sectors 
around the points of (, and it can be computed with an arbitrarily high precision. 
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• In the case of a Siegel disc D, information that would allow us to identify a repelling 
orbit Cd in the same connected component of Kp as D. Such an orbit always exists, 
and can be identified using a finite amount of combinatorial information. 

2.1. Computing Kp. We are given a dyadic point rf G © and an n G N. Our goal is to 
always terminate and output 1 ii B{d,2~'^)r\Kp ^ and to output ii B{d,2-2'^)nKp = 0. 
We do it by constructing five machines. They are guaranteed to terminate each on a different 
condition, always with a valid answer. Together they cover all the possible cases. 

Lemma 2.1. There are five oracle machines M^xt, Mjui, Mhyp, Mpar, Mgieg such that 

(1) if d is at distance > | ■ 2~" from Kp, Mf.xt{d,n) will halt and output 0. If d is at 
distance < 2"" from Kp, M^xtidi n) will never halt; 

(2) if d is at distance < | ■ 2^" from Jp, Mjui{d,n) will halt and output 1. If d is at 
distance > 2 ■ 2"" from Jp, Mjui{d,n) will never halt; 

(3) Mfiyp{d,n) halts and outputs 1 if and only if d is inside an attracting basin of a 
hyperbolic orbit of p; 

(4) Mpar{d,n) halts and outputs 1 if and only if d is inside an attracting basin of a 
parabolic orbit of p; 

(5) Msieg{d,n) halts and outputs 1 if the orbit of d reaches a Siegel disc, and d is at 
distance > | ■ 2~" from Jp. It never halts if d is at distance > 2 ■ 2"" from Kp. 

Recall the Fatou-SuUivan classification of Fatou components of a polynomial mapping of C 
(see e.g. JMlIl): 

Theorem 2.2 (Fatou-Sullivan classification). Every Fatou component of a polynomial 
mapping of C of degree at least two is a preimage of a periodic component. Every periodic 
component is of one of the following types: the immediate basin of an attracting (or a super- 
attracting) periodic point; a component of the immediate basin of a parabolic periodic point; 
a Siegel disk. 

Proof of Theorem \1.4[ given Lemma \2.1\ By Fatou-Sullivan classification it is not hard to 
see that for each (d, n) at least one of the machines halts. Moreover, by the definition of 
the machines, they always output a valid answer whenever they halt. Hence running the 
machines in parallel and returning the output of the first machine to halt gives the algorithm 
for computing Kp. D 

It remains to prove Lemma f2. II 

Proof, (of Lemma 12. 1|) We give a simple construction for each of the five machines. 

(1) Mea;^: Take a large ball B such that Kp C B. Intuitively, we pull the ball back under 
p to get a good approximation of Kp. Let B^ be a 2~'^"+^)-approximation of the set 
p'^{B). Output iff -Bfc n -B(ci, I ■ 2"") = 0. It is not hard to see that this algorithm 
satisfies the conditions on M^xt- 

(2) MjuV- Enumerate all the repelling periodic orbits of p. Let C^ be a 2~'^"+^)-approximation 
of the union of the first k orbits enumerated. Output 1 iff d{d, Ck) < ^ ■ 2"". The 
repelling periodic orbits are all in Jp and are dense in this set. Hence the algorithm 
satisfies the conditions on Mjui- 
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(3) Mhyp'. Let Zk be a 2~'^-approximation oi p''{d). If d isinside the basin of attraction 
for a some orbit (, then Zk for some k will be inside Z)^ for some k. Output 1 if Zk is 
at least 2~^-far from the boundary of D^. 

(4) Mpar'- Very similar to M^yp. The only difference is that now we are checking for 
convergence to an attracting petal of a parabolic orbit. 

(5) Msi(.g\ This is the most interesting case. It is not hard to see that for each fc, we can 
compute a union Ek of dyadic balls such that 

[jp^ (B{d, ^ ■ 2-")] CE,c[jp' (sid, ^ ■ 2-")] . 
j=0 ^ ^ j=0 ^ ^ 

Let c be the center of the Siegel disc (one of the centers, in case of an orbit), and let 
y be the given periodic point in the connected component of c. We terminate and 
output 1 if En separates c from y in C (or covers either one of them) . 

If d is inside the Siegel disc, then the forward images of B{d,^ ■ 2^") will cover 
an annulus in the disc that will separate c from the boundary of the disc, and in 
particular from y. Hence Mgieg will terminate and output 1. 

If the distance from d to K^ is > 2 ■ 2~", then Ej. fl Kp = for all k. In particular, 
Ek cannot separate c from y, since they are connected in Kp. 

D 

2.2. The quadratic case. In the quadratic case there is at most one non-repelling orbit. 
In the case there is a Siegel disc, the Julia set is connected, and any repelling periodic orbit 
can be taken as the orbit connected to the Siegel disc. In fact, it is not hard to see that if we 
exclude the parabolic case, one machine suffices to take care of all the connected filled Julia 
sets. As a corollary we get: 

Corollary 2.3. Denote by M the Mandelbrot set, and by P the set of c's for which Jc is 
parabolic. The function K : c ^-^ K^^^c is continuous in the Hausdorff metric on the set 
M -P. 

3. COMPUTABILITY OF THE SET J 

Recall that 

J = {(z,c) : zeJc} C C^. 
Theorem 11.51 asserts that J is computable. We prove it by showing that J is weakly com- 
putable. 

Definition 3.1. We say that a set C is weakly computable if there is an oracle Turing 
Machine M'^{n) such that if represents a real number x, then the output of M'^{n) is 

( 1 iixeC 

M't'{n) = l if S(x, 2-") n C = 

I or 1 otherwise 

It has been shown that the weak definition is equivalent to the standard definition. See 
|Brv2j ■ for example. We will need the following lemma. 

Lemma 3.1. For any point {z,c) in the complement of the closure J, z converges to an 
attracting periodic orbit of fc '. z ^—>- z"^ + c. 



The proof of the lemma occupies 21 

The following lemma allows us to "cover" all points that belong to J. 

Lemma 3.2. There is an algorithm Ai{n) that on input n outputs a sequence of dyadic 
points pi,p2, ■ ■ ■ & C X C such that 

oo 

i=i 

Proof. It is well known that repelling periodic orbits of fc are dense in Jc- Hence, the set 

Srep = {(^? c) : z is in a repelling periodic orbit of fc} 

is dense in J. Srep is a union of a countable number of algebraic curves 5*^^ given by the 
constraints 

/rw = ^ 
i(/r)'(^)i>i 

For each m we can compute a finite number of points p"^, ■ ■ ■ ,p5?^ approximating S^p such 
that 

S(5,^p,2-("+=^)) C |Jfi(pf,2-("+2)) c S(5,":p,2-("+^)). 
i=i 
We have 

oo 

J Orep II i^rep' 

m=l 

Hence the computable sequence p}, . . . ,pj^,p^, . . . ,^^^5 • • • ^P™) • • • ;Prvi' • • • satisfies the con- 
ditions of the lemma. D 

Corollary 3.3. There is an oracle machine M^ (n), where (pi is an oracle for z G C and 
02 is an oracle for c G C, such that Mf'^''''^ always halts whenever d{{z,c),^) < 2^^"+^^ and 
never halts ifd{{z,c),S) > 2^". 

Proof. Query the oracles for a point p E C x C such that d{p, {z,c)) < 2^^"'+'^^ Then run 
the following loop: 

i^O 
do 

i ^i + 1 

generate pi using Ai{n) from Lemma f3. 21 

while d{p,pi) > 2-("+2) 

If d{{z,c),S) < 2"("+^\ then d{p,S) < 2"^'^+^^ hence by Lemma IT^ there is an i such 
that d{p,pj) < 2~("+^), and the loop terminates. If (i((z,c), J) > 2"", then (i(p, J) > 2"" - 
2~(n-4) > ^5 . 2-(n+i). Hence, by Lemma O P i fife, 2-("+i)) for all i, and the loop will 
never terminate. D 

The following lemma allows us to exclude points outside J from JI. 

Lemma 3.4. There is an oracle machine M2 , where (pi is an oracle for z E C and (p2 is 
an oracle for c E C, such that M2 halts if and only if z converges to an attracting periodic 
orbit (or to 00^ under fc'.z^^z^ + c. 



Proof. M2 is systematically looking for an attracting cycle of fc- It also iterates fc on z with 
increasing precision and for increasingly many steps until we are sure that either one of the 
two things holds: 

(1) the orbit of z converges to 00; or 

(2) we find an attracting orbit of fc and the orbit of z converges to it. 

If the search is done systematically, the machine will eventually halt if one of the possibil- 
ities above holds. It obviously won't halt if neither holds. D 

Proof, (of Theorem II. 5j) The algorithm is: Run the machines Mf^ (n) from Corollary 
13.31 and M2 from Lemma 13.41 in parallel. Output 1 if Mi terminates first and if M2 
terminates first. 

First we observe that Mi(n) only halts on points that are 2~"-close to J, in which case 
1 is a valid answer according to Definition 13.11 Similarly, M2 only halts on points that are 
outside J, in which case is a valid answer. Hence if the algorithm terminates, it outputs a 
valid answer. It remains to see that it does always terminate. Consider two cases. 

Case 1: (^, c) G J. In this case d{{z,c),3i) = < 2~^"+^\ and the first machine is 
guaranteed to halt. 

Case 2: {z,c) ^ J. By Lemma l3.H z converges to an attracting periodic orbit of fc in 
this case, and hence the second machine is guaranteed to halt. D 

4. Proof of Lemma IXTl 

Suppose z ^ Jc and the orbit of z does not belong to an attracting basin. By the Fatou- 
Sullivan classification (see e.g. |Milj ) . there exists A; G N such that w = fc{z) belongs to a 
Siegel disk or to the immediate basin of a parabolic orbit. Our aim is to show that for an 
arbitrary small S > 0, there exists a pair (z, c) G C x C with \z — z\ < 6, \c — c\ < 6, and for 
which z G Jc- We will treat the Siegel case first. 

4.1. The case when w lies in a Siegel disk. Let us denote A the Siegel disk containing 

w, and let m G N be its period, that is, the mapping 

/r : A ^ A 

is conjugated by a conformal change of coordinates : A ^ D to an irrational rotation of 

D. 

The following statement is elementary (cf. Prop. 7.1 in |Douj ) : 

Proposition 4.1. Denote ( = 0^^(O) G A the center of the Siegel disk. For each s > 
there exists c G B{c, s) such that fc has a parabolic periodic point ( of period m in B{(, s). 
In particular, Jc is connected, and B{(, s) Cl Jc ^ 0. 

Consider now the /"-invariant analytic circle 

^^ = 0-i({;2 = re2-^, ee [0,27r)}) 
which contains w. Let e > be such that 

B{w,e)cf^{B{z,6))nA. 



Set B = B{w, e/2) and let n G N be such that the union 

0<i<n 

By Proposition l4. ll for all (5 > small enough, there exist c G B{c, S) for which Jg is connected 
and there is a point of Jg inside the domain bounded by Sr- Since repelling periodic orbits 
of fc are dense in dA, again for 6 small enough, there are points of Jg on the outside of Sr as 
well, and so there exists a point ^ G Jg fl S",.. By construction, there exists j G N such that 
f^{B{z, 6)) 3 C,. By invariance of Juha set, if c is close enough to c we have B{z, 5)\~\Jc^ 0, 
and the proof is complete. 

4.2. The case when w lies in a parabolic basin. We need to recall the Douady-Lavaurs 
theory of parabolic implosion f |Dou|lLavj ). Denote C the parabolic periodic point of fc whose 
immediate basin contains w, and let m G N be its period. 

Recall that an attracting petal Pa is a topological disk whose boundary contains C, such 
that fl^^{PA) C Pa for some /c G N, and such that the quotient Riemann surface 

Ca = PaIIT^ ^ c/z. 

The quotient Ca-, called an attracting Fatou cylinder, parametrizes the orbits converg- 
ing under the dynamics of the iterate f^ to C inside the periodic cycle of petals Pa, 
/^(Pa), ■ ■ ■ , fc (Pa)- Recall (see |Milj ) that a quadratic polynomial fc has only one 

cycle of petals at (. A repelling petal Pr is an attracting petal for the local inverse /"*" 
fixing (; the union 

U friPAUPn) 

0<i<k-l 

forms a neighborhood of (. The repelling Fatou cylinder Cji is defined in a similar fashion. 
Let r be any conformal isomorphism Ca -^ Cr. After uniformization, 

Ca ^ C/Z, Cr ^ C/Z 

t{z) = z + q mod Z for some q E C Let g : Pa ^ Pr he any lift of r; it necessarily 
commutes with /^^. Consider the semigroup G generated by the dynamics of the pair 
{fc,g)- The orbit Gz of a point 2; G C is independent of the choice of the lift g and only 
depends on r. 

Set 

■Jic,T) = {z EC such that GzHJc^^}- 
It can be shown that this set is the boundary of 

K{c,t) = {z E C such that Gz is bounded}. 

The Douady-Lavaurs theory postulates: 

Theorem 4.2. For every r as above and every s > there exists c E B{c,s) such that 
B{Jc,s) D J^c,t)- 

Since ( E Jc, and Jc is connected, there exists a point m G Jc fl Pr. Let w E Ca be the 
orbit of w, and let u E Cr be the orbit of u. Choose r : Ca — ^ Cr so that t{w) = u. Then 
J{c,t) 3 -2, and the claim follows by Theorem 14.21 
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